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Abstract 

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases 
affecting humans and animals. The neurodegenerative diseases such as Creutzfeldt- 
Jakob disease, variant Creutzfeldt-Jakob diseases, Gerstmann-Straussler-Scheinker 
syndrome, Fatal Familial Insomnia, Kuru in humans, scrapie in sheep, bovine spongi- 
form encephalopathy (or mad-cow disease) and chronic wasting disease in cattle etc 
belong to prion diseases. By now there have not been some effective therapeutic ap- 
proaches to neat all these prion diseases. In 2008, canine mammals including dogs 
(canis familials) were the first time academically reported to be resistant to prion dis- 
eases (Vaccine 26: 2601-2614 (2008)). Rabbits are the mammalian species known to 
be resistant to infection from prion diseases from other species (Journal of Virology 77: 
2003-2009 (2003)). In 2010 horses were reported to be resistant to prion diseases too 
(Proceedings of the National Academy of Sciences USA 107: 19808-19813 (2010)). 
By now all the NMR structures of dog, rabbit and horse prion proteins had been re- 
leased into protein data bank respectively in 2005, 2007 and 2010 (Proceedings of the 
National Academy of Sciences USA 102: 640-645 (2005), Journal of Biomolecular 
NMR 38:181 (2007), Journal of Molecular Biology 400: 121-128 (2010)). Thus, at 
this moment it is very worth studying the NMR molecular structures of horse, dog and 
rabbit prion proteins to obtain insights into their immunity prion diseases. 

This article reports the findings of the molecular structural dynamics of wild-type 
horse, dog, and rabbit prion proteins. The dog and horse prion proteins have stable 
molecular structures whether under neuttal or low pH environments. Rabbit prion pro- 
tein has been found having stable molecular structures under neutral pH environment, 
but without sttuctural stability under low pH environment. Under low pH environment, 
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the salt bridges such as D177-R163 were broken and caused the collapse of the stable 
a-helical molecular structures. The contribution of D177-R163 salt bridge was con- 
firmed by the recent experimental results on the X-ray and NMR structures of rabbit 
prion protein (Proceedings of National Academy of Sciences USA 107: 19808-19813 
(2010), Journal of Biological Chemistry 285 (41), 31682-31693 (2010), Journal of 
Molecular Biology 400: 121-128 (2010)). This salt bridge might have a value to the 
scientific community in its drive to find treatments for prion diseases. 

Keywords: Prion diseases; Immunity; Rabbit, dog, and horse prion proteins; Molecular 
dynamics. 



1. Introduction 

Prion diseases, including Creutzfeldt-Jakob disease (CJD), variant Creutzfeldt-Jakob dis- 
eases (vCJD), Gerstmann-Straussler-Scheinker syndrome (GSS), Fatal Familial Insom- 
nia (FFI), Kuru in humans, scrapie in sheep, bovine spongiform encephalopathy (BSE 
or 'mad-cow' disease) and chronic wasting disease (CWD) in cattle, etc., are invari- 
ably fatal and highly infectious neurodegenerative diseases affecting humans and ani- 
mals. However, for treating all these diseases, there is no effective therapeutic approach 
| Aguzzi et al. 2006[lPrusiner 1998UWeissmann 20041. 



Fortunately, dogs, rabbits, and horses have been reported to be the few mammalian 
species being resistant to infection from prion diseases. Their NMR molecular structures 



were released into the Protein Data Bank [Berman et al., 2000| in the years 2003, 2005, 



2010 respectively. The N-terminal residues of prion proteins are unstructured and not in- 
cluded in this study. The C-terminal residues are well structured, with 3 a-helices, 2 short 
anti-parallel /3-sheets, and a disulfide bond between the 2nd and 3rd a-helix (residues num- 
ber 178 and number 213). The infectious prion (PrP Sc ) is an abnormally folded form of 
the normal cellular prion (PrP c ) and the conversion of PrP c to PrP Sc is believed to involve 
conformational change from a predominantly a-helical protein (42% a-helix, 3% /3-sheet) 
to a protein rich in /3-sheets (30% a-helix, 43% /3-sheet). 

Rabbits are one of the few mammalian species reported to be resistant to infection 



from prion diseases isolated from other species [Vorberg et al., 2003]. Recently, the NMR 
molecular structures of wild-type, mutant S173N and mutant 12 14V rabbit prion proteins 
(124-228) were released into the Protein Data Bank |Berman et al, 2000] with PDB ID 
codes 2FJ3, 2JOH, 2JOM respectively. Zhang (2010a) studied these NMR structures by 
MD simulations and simulation results at 450 K under low and neutral pH environments 
confirmed the structural stability of wild-type rabbit prion protein. However, the previous 
MD simulation results at 500 K under neutral pH environment showed that wild-type rabbit 
prion protein (124-228) does not have a structural stability [Zhang 2009 ]. At 450 K, Zhang 
(2011a) also studied the molecular dynamics of NMR structures of human and mouse prion 
proteins (lQLX.pdb and lAG2.pdb respectively) to make comparisons with that of rabbit 
prion protein [Zhang 2011a]. The comparisons show that the enhanced stability of the C- 
terminal ordered region especially helix 2 through the D177-R163 salt-bridge formation 
renders the rabbit prion protein stable, the salt bridge D201-R155 linking helices 3 and 1 
and the hydrogen bond H186-R155 also contribute to the structural stability of rabbit prion 
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protein; however, the human and mouse prion protein structures were not affected by the 
removing these two salt bridges. 

Dogs (canis familials, canine mammals) were the first time academically reported in 
2008 to be resistant to prion diseases [Polymenid oua et al. 2008[ . Zhang et al. (201 1) stud- 
ied the molecular dynamics of dog prion protein ( lXYK.pdb) at 450 K and found that the 
dog prion protein has stable molecular structures whether under neutral or low pH environ- 
ments [Zhang 201 lb | . 

Horses were academically reported to be resistant to prion diseases recently 
UPerez et al. 20111 . Zhang (2011) studied the molecular dynamics of horse prion protein 
(2KU4.pdb) at 450 K and found that the dog prion protein has stable molecular structures 
whether under neutral or low pH environments [Zhang 201 lc[. Molecular dynamics results 
show that the strong salt bridge ASP177-ARG163 contributes to the structural stability of 
horse prion protein [Zh ang 2011c| . 

All the above molecular dynamics simulations done for rabbit, dog and horse prion 
proteins are at 450 K and 300 K under neutral and low pH environments. The simu- 
lations done at room temperature 300 K whether under neutral or low pH environment 
display very little fluctuation for rabbit, dog, horse, human and mouse prion proteins. 
The common result for rabbit, dog and horse prion proteins found at 450 K is: there 
is a strong bridge ASP177-ARG163 (like a taut bow string) keeping the (32-a2 loop of 
prion proteins linked. The finding agrees with the recent experimental results on the j32- 
a2 loop IBett et al, 2012] ILin et al. 20m IFernandez-Funez et al. 20TT1 IPerez et al. 20TT1 
IKhan et al. 20T0l |Sweeting et al. 2009} IWen et al. 20T0al I Wen et al. 20T0bl . which plays 
an important role to stabilize the structural stability of prion proteins. 

This Chapter reports the molecular dynamics simulation results done for rabbit, dog 
and horse prion proteins are at 350 K under neutral and low pH environments. Molecular 
dynamics results agree with the results found at 450 K. The stable a-helical secondary 
structure of rabbit prion protein collapses under low pH environment. The strong bridge 
ASP177-ARG163 (like a taut bow string) keeping the /32-a2 loop of prion proteins linked 
contributes to the structural stability of rabbit, dog and horse prion proteins. This salt bridge 
might be a potential drug target and have a value to the scientific community in its drive 
to find treatments for prion diseases. The rest of this article is arranged as follows. The 
molecular dynamics simulation materials and methods are introduced in Section 2. Section 
3 mainly gives molecular dynamics simulation results and their discussions. Concluding 
remarks are given in the last Section. 



2. Materials and Methods 

The molecular dynamics simulation materials and methods for horse, dog and rabbit prion 



proteins are completely same as the ones of [Zhang 2010 1. Simulation initial structure 
for the dog and rabbit prion proteins were built on HoPrP c (l 19-231) (PDB entry 2KU4), 
DoPrP c (121-231) (PDB entry 1XYK) and RaPrP c (124-228) (PDB entry 2FJ3), respec- 
tively. Simulations were done under low pH and normal pH environments respectively. 
All the simulations were performed with the AMBER 11 package BCase et al. 20101 . with 
analysis carried out using functionalities in AMBER 9 ICase et al. 20061 and AMBER 7 
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CARNAL llCase et al. 20021 . Graphs were drawn by XMGRACE of Grace 5.1.21, DSSP 
Kabsch et al. 19831 . 

All simulations used the ff03 force field of the AMBER 11 package, in neutral and 
low pH environments (where residues HIS, ASP, GLU were changed into HIP, ASH, GLH 
respectively by the XLEaP module of AMBER 11 in order to get the low pH environment). 
The systems were surrounded with a 12 angstrom layer of TIP3PBOX water molecules and 
neutralized by sodium ions using XLEaP module of AMBER 1 1 . The solvated proteins 
with their counterions were minimized mainly by the steepest descent method and then a 
small number of conjugate gradient steps were performed on the data, in order to remove 
bad hydrogen bond contacts. Then the solvated proteins were heated from 1 00 to 300 
K during 1 ns (with step size 1 fs) and from 300 to 350 K during 1 ns (with step size 
2 fs). The thermostat algorithm used is the Langevin thermostat algorithm in constant 
NVT ensembles. The SHAKE algorithm and PMEMD algorithm with nonbonded cutoffs 
of 12 angstrom were used during the heating. Equilibrations were done in constant NPT 
ensembles under Langevin thermostat for 2 ns. After equilibrations, production MD phase 
was carried out at 350 K for 30 ns using constant pressure and temperature ensemble and 
the PMEMD algorithm with nonbonded cutoffs of 12 angstrom during simulations. Step 
size for the production runs is 2 fs. The structures were saved to file every 1000 steps. 



3. Results and Discussion 



At 350 K there are fluctuation and variation among rabbit, dog and horse prion proteins, but 
their backbone atom RMSDs (root mean square deviations) respectively calculated from 
their minimized structures and their radii of gyrations do not have great difference even 
under low pH environment. Their secondary structures under neutral pH environment at 
350 K do not change very much either. However, the secondary structures under low pH 
environment at 350 K have great differences between rabbit prion protein and dog and horse 
prion proteins (Figure [T]): the a4ielices of rabbit prion protein were completely unfolded 
and began to turn into /3-sheets but those of dog and horse prion proteins were not changed 
very much. These results indicate the C-terminal region of RaPrP c has lower thermosta- 
bility than that of DoPrP c and HoPrP c . Under the low pH environment, the important salt 
bridges such as D177-R163, D201-R155 were removed (thus the free energies of the salt 
bridges changed the thermostability) so that the structure nearby the central helices 1-3 was 
changed for rabbit prion protein. 

There always exists a strong salt bridge (where the oxygen-nitrogen distance cut-off 
calculated for the salt bridges is 3.2 angstroms) between D177-R163 for rabbit prion pro- 
tein, and D178-R164 for dog and horse prion proteins (Figure [2]) under neutral pH en- 
vironment at 300 K, 350 K, and 450 K. We can see in Figure [2] that this salt bridge is 
just like a taut bow string linking /32-sheet and a2-helix. It contributes to the struc- 
tural stability of all these prion proteins. Thus, the region of (32-a2 loop should be 
a potential drug target region according to recent reports on the /32-a2 loop of prion 
proteins IBett et al, 2 012J ILin et al. 20111 IFernandez-Funez et al. 20111 IPerez et al. 20111 
IKhan et al. 20T0l |Sweeting et al. 2009[IWen et al. 20"T0allWen et al. 2010bH . 

Lastly, we compare rabbit with dog and horse prion proteins through the pairwise se- 
quence and structure alignments. We find that there are more identities and similarities be- 
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Figure 1. Secondary structures of rabbit, dog and horse prion proteins (from up to down) at 
350 K under neutral to low pH environments (from left to right) (X-axis: ns - 30 ns (from 
left to right), Y-axis: residue numbers 124 - 228 / 121 - 231 / 119 - 231 (from up to down)). 




Figure 2. The salt bridge D177-R163 of rabbit prion protein, D178-R164 of dog and horse 
prion proteins. 
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tween horse and dog prion proteins, compared with the identities and similarities between 
dog and rabbit prion proteins, and between horse and rabbit prion proteins. This point can be 
furthermore confirmed by the following data of the Needleman-Wunsch or Smith- Waterman 
Pairwise Sequence Alignment and the jCE algorithm or jCE Circular Permutation Pairwise 
Structure Alignment. (1) For dogs and horses, the Identities are 89.38 % (query) and 90.99 
% (subject) and the Similars are 93.81 % (query) and 95.50 % (subject) for the Sequence 
Alignment, and the Identity is 89.1 1% and the Similarity is 93.07% for the Structure Align- 
ment. (2) For dogs and rabbits, the Identities are 90.09 % (query) and 72.46 % (subject) 
and the Similars are 95.95 % (query) and 77.18 % (subject) for the Sequence Alignment, 
and the Identity is 92.93 % and the Similarity is 98.99 % for the Structure Alignment. (3) 
For horses and rabbits, the Identities are 87.61 % (query) and 71.74 % (subject) and the 
Similars are 93.81 % (query) and 76.81 % (subject) for the Sequence Alignment, and the 
Identity is 71.72% and the Similarity is 77.78% for the Structure Alignment. This might 
explain the reasons why rabbit prion protein differs very much from dog and horse prion 
proteins in secondary structures under low pH environment [Zhang 201 lb, Zhang 2 01 lc[ . 



4. Conclusion 

The paper is a straight forward molecular dynamics simulation study of rabbit, dog and 
horse prion proteins which apparently resist the formation of the scrapie form. The analyses 
of molecular dynamics results confirmed the structural stability of rabbit prion protein under 
neutral pH environment, and of dog and horse prion proteins under neutral and low pH 
environments. The main point is that the enhanced stability of the C-terminal ordered region 
especially helix 2 through the D177-R163 (D178-R164) salt-bridge formation renders the 
rabbit, dog and horse prion proteins stable. This salt bridge might be a potential drug target 



for prion diseases. Zhang (2011) systemically analyzes this salt bridge in [Zhang 201 Id | 
and, like GN8 of [ [Kuwata et al., 20071 , fixing the distance between D177 and R 163 to form 
a salt bridge distance might be able to design a drug for treating prion diseases. 
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